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Abstract 



We consider a generalisation of the stochastic Radon transform, introduced for an in- 
verse problem in tomography by Panaretos. Specifically, we allow the distribution of the 
three-dimensional rotation in the statistical model of that work to be different from Haar 
measure, and to possess only the weaker symmetry property of conjugation-invariance. As 
a preparation we study, based on simple geometric ideas, the relationship between moments 
of two univariate densities describing conjugation-invariant rotations. For the class of ro- 
tations introduced by Leon, Masse and Rivest, the generalised stochastic Radon transform 
turns out to have a 'fake uniformity' property. We also derive the binary Bayes classifier 
l/*~) for rotations whose distributions are deterministic shifts of a common conjugation- invariant 

(N 



distribution. 
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1 Introduction 

Random rotations in Euclidean space R (d > 2) arise in a variety of contexts in stochastic geom- 
etry, especially in the construction of other stochastic models, for example random hyperplane 
processes; see the monograph [13]. In these contexts, the rotation distribution is usually taken as 
the Haar probability measure \i = \i& on the d-dimensional rotation group SO(d). Recall that jjl 
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is characterised by the property that fi(QA) = fi(AQ) = fi(A) for any Q e SO(d) and any mea- 
surable A C SO(d), where 50(d) is endowed with its Borel cr-field, QA = {QP : P G SO(d)}, 
and AQ is defined analogously. See [13, Ch. 13] for more information and one of several con- 
structions of Haar measure on SO(d). Most applications concern the case d = 3, not only in 
stochastic geometry but also statistics. Two important instances are a statistical model for 
random orientations of limbs of experimental subjects performing drilling tasks, studied in [4], 
and a model for tomographic imaging introduced in [7] and further studied in [8]. In [4], a 
class of rotation distributions in R was introduced that was shown to be amenable to exact 
computations of several key quantities including first and second moments. In [7], the goal 
was to estimate the shape (as defined in that paper) of a particle imaged with cryo-electron 
microscopy, where the orientation of the particle is not only random but also unknown, its 
orientation being governed by /x. As was shown in that paper, this last assumption is necessary 
in order to produce a well-defined estimator of the shape of the particle. 

The basis of the present paper is a synthesis of [4] and [7]. Our focus is on probabilistic rather 
than statistical aspects of the model. Specifically, we study a probability model that arises 
when the rotation distribution in the stochastic Radon transform introduced in [7] is not Haar 
measure, but only possesses the more general symmetry property of conjugation-invariance, 
whose definition we recall in Section 2. Our results will be more general than what we need 
for the tomography application, and in particular, we will obtain several new results relating 
to three-dimensional rotation distributions. We also derive the Bayes classification rule with 
two classes which are shifts of a common conjugation-invariant distribution, and examine the 
behaviour for small shifts. For the case of two multivariate Normal distributions on M d with 
equal covariance matrices, a corresponding result is well known; see [10, pp. 21-22]. As far as we 
are aware, our approach is essentially new, although there is some similarity in some instances 
with the geometric reasoning in [3], which studied the distribution of the projected shape under 
Haar measure. 

Our two main results, namely Theorems 1 and 2 in Sections 4 and 5 respectively, yield two 
different ways of quantifying the effect of moving away from the Haar measure situation. In 
relation to the first of these results, a finding of particular interest is what we call a fake 
uniformity property, 'uniform' being taken here and below as a synonym for Haar measure \x 
(see the discussion around Definition 1 in Section 4). The second result quantifies how well one 
may distinguish two distributions for a shift of 'modal' angle, as explained later; as expected, 
we find better accuracy for stronger deviation of the common 'centred' distribution from \jl 
(see the discussion around (12) in Section 5). We shall not attempt to develop a variational 
approach using distances between rotation distributions as was done in [9], as this seems a 
secondary matter to the phenomena we study. Our results also provide a new interpretation of 
the concentration parameter in statistical models for random rotations. Finally, our results may 
provide a foundation for statistical tests that are based on the quantities we consider. It seems 
clear that the fake uniformity property is a problem for the most basic of such tests, namely 



that for uniformity of the rotation. 

Section 2 contains background on rotations, as well as notation. Section 3 investigates the 
connection between two univariate distributions which describe conjugation- invariant rotations. 
Section 4 contains the main result on the stochastic Radon transform for what we call shift- 
symmetric rotations, along with an example that reveals the fake uniformity property and 
motivates us to define it separately. Finally, in Section 5 the Bayes classifier is derived. 



2 Notation and Background 

All random variables and random elements on other state spaces are taken to reside on a common 
probability space (fi, & , P), with E(-) denoting the corresponding expectation operator. 

We shall be concerned not only with random variables on M, but also on the sphere in three 
dimensions § 2 := {i £ I 3 : ||x|| = 1}, and on 50(3). If T denotes a random element in either 
of these three spaces, we shall denote its density by /t, it being understood that the dominating 
measure for this density is ordinary Lebesgue measure, normalised Lebesgue surface measure a 
on S 2 , or Haar probability measure fi on 5*0(3) respectively. As done above, we use boldfacing 
to distinguish deterministic and random rotations, but will use uppercase and lowercase letters 
for vectors and scalars as is customary. All vectors are regarded as column vectors. We write 
Vi = v ■ e^> for v £ M 3 , where the dot denotes the usual inner product and e^' is the ith column 
of the 3x3 identity matrix, for i = 1,2,3. The matrix diag(ai, 02, 03) is the 3x3 diagonal 
matrix with entries a\, 02, 03 along its main diagonal. The Beta integral is 



B(p,q)= f xP-\l-xf-\ P,q>0. 

Jo 



and Beta(p, q) is the Beta distribution with parameters p, q > 0. Recall that if W ~ Beta(p, q) 

then 

E(W r ) = ^ r for any integer r > 0, (1) 

[p + q)r 

where (x) n = r(x + l)/r(x — n + 1) is the Pochhammer symbol and T(-) is the Gamma function. 
We use the notation E(A; B) := E(Xl^) for a scalar random variable X and a measurable set 
Bel, where 1b is the indicator function of B. 

Recall that any rotation P G 50(3) (except if it belongs to the null set comprising the identity 
and rotations by an angle n) may be expressed uniquely in terms of the rotation axis u € S 2 



and rotation angle 9 G [0,7r), as follows (see [11, pp. 188-189] or [12, p. 99]): 

/o -, 

p = cos(0)I + sin(0)<S(u) + (1 - cos(6))uu T , where 5(a) = a 3 -oi I • (2) 




We assume throughout that the letter R denotes a random three-dimensional rotation which 
is conjugation-invariant, meaning that is satisfies R = Q T RQ, where = denotes equality in 
distribution; then automatically R = R. The terminology here follows [12]. We say that R 
(or its distribution) belongs to Z s . This definition for Z s is somewhat less restrictive than the 
one given in [11, p. 190]. In the context of random P at (2), one has fp G X s if and only if 
U ~ 02 (cf. [11, p. 191]). In the sequel, when referring to random elements U and 0, we always 
assume them to be defined through P via the bijection (up to a null set) P -H- (U, ©). 

To date, only few rotation distributions have been studied in the literature. We shall consider 
the distribution introduced in [4], referred to below as the Cayley-LMR distribution: 

/CAY(M, K )(P) = ^r{n + iri) {l + iT{PMT)) ^ (3) 

and the symmetric matrix Fisher-von Mises distribution, whose density is proportional to 
ex.p{n(PM T )}; see [5, pp. 289ff] for more information on the latter distribution. In both cases, 
k > is a dispersion parameter, with n = corresponding to the case fp = fn = /x, and M 
may be regarded as a modal matrix. We shall make the minimal assumption that our rotation 
distributions are parametrised by modal matrices in the sense of a group model, see [9, p. 653], 
and a dispersion parameter k, with the case k = corresponding to Haar measure [i. They are 
allowed to have other parameters as well. 

We shall on several occasions consider the random variable 

1 + cos e ( 2 ( 6 



and its (central) moments {p r = E(X r ), r > 0}. For the Cayley-LMR distribution, X ~ 
Beta(/s + 1/2, 3/2) ([4], p. 425) and thus by (1), 

k + 1/2 (« + l/2)(« + 3/2) 

Pl(K) = TT^' P2{K)= ( K + 2)(k + 3) • (5) 

It will become clear later why we regard these as functions of k. 

Finally, we introduce some convenient terminology. For a random rotation P which is such that 
there is a conjugation-invariant rotation R G SO (3) and a deterministic rotation M £ SO(3) 
such that PM T = R, we write P G X^ . We call Q shift-symmetric, and write Q G Z^, if there 
exists an M such that Q G Z^. The same terminology applies to the underlying distributions 



of the random variables. 



3 Relation between two distributions and moments 



In this section we investigate the connection between two key distributions derived from a 
three-dimensional rotation. The first proposition shows that if fa G X s then the density /r 6 (3) 
is simply that of 

Z=(e^fRe^. (6) 

In the proof and in several later instances we shall need the following simple result on integration 
on the sphere, which is a special case of [1, Proposition 9.1.2, p. 189]: 

Lemma 1. Let / be an integrable function on § 2 which is zonal, that is, f(x) = F(xs), where 
F is a function defined on [—1,1]. Then 

1 '- 1 



f(x)a(dx) = - / F(t)dt. 

In particular, if U ~ oi then U% ~ Uniformf— 1, 1]. 

Proposition 1. fg £ (z){w) = fziwz)- Moreover, / Be (3) determines /r uniquely. 

Proof. Consider realisations R = R(ojq), z = Z(ujq), u = U(ljq), and 9 = Q(ojo). By the law of 
cosines for spherical triangles (see [2, p. 54]), we have 



z = cos 



:(e (3) ,i?e (3) )) = cos 2 (Z(u,.Re (3) ))+sin 2 (z(u, Re {3) )) cos^ = u§+(l-«§)cose, (7) 



where Z.(v,w) = arccos(f • w) G [0, 7r) is the angle between the vectors v,w G S 2 . Hence 
cos0 = (Z — f7 2 )/(l — C/|) almost surely. The equality of the densities follows from Lemma 1, 
cf. also the remark in [1, p. 189]. Now (2) shows that the distribution of R, and thus also that 
of Re^ 3 ' , is determined uniquely. □ 

Corollary 1. If Q € Xf with M such that Me' 3 ' = e( 3 ), then f Qe3 (w 3 ) = fz{w 3 ). 

Proof. Follows from the fact that [9, p. 653] fQ ea = file?, ° n _1 f° r an Y n G 50(3) with 

Me( 3 '=n e ( 3 ). □ 



Proposition 2. 



f (s) = — I ^ M=L dx -1< 8 < 1 

JZ[S) 2^2 Jo y / (l + S-2x)(l-x) X} " S ^ 



Proof. By Lemma 1, (7) and Fubini's theorem, 

P(Z<s) = - J P(Z < s | U 3 = t)dt= - I P{t 2 + (l-t 2 )cos9 < s}dt 
2 J -i 2 J_ 1 

l+s 



I f^r ll + 8-2x 

fx{x)dx, 



y/2 J V 1 - x 

and differentiation with respect to s completes the proof. □ 

Remark 1. Conversely, it is readily shown that 

i r 1 

fx{s) = - f z {2s-2st + 2t 2 -l)dt, 0<s<l. 

The relationship between the densities in Proposition 2 will now be used to relate the corre- 
sponding moments up to some finite order. Consider the vector r = (tj)' 1 ?^, m > 1, where 
Tj = E(/P). Lemma 1 and Proposition 2 yield for k > 1, by partial integration, 

n _ J_ /' t * dt /* /*(*> dl 

2\/2 y.i A N /(l + (-2l)(l-l0 

1 rhto.tef 1 j a, 



2^/2 7o Vl-x y 2 *-i \/l + t-2x 

= 1 _kf 1 Mx^ dx r 1 t k-i Vl + t _ 2xdt 

V2 Jo V 1 - X ./2x-l 

Another partial integration gives, for any /c > 1 and — 1 < x < 1, 
G k (x) := /"' t fe -Vl + t-2xdt = ^(1 - x) 3 / 2 - fc^ /"' t k - 2 (l + i - 2x) 3 / 2 dt 

J2x-1 3 3 72x-l 

= afc(x) + fe fe (x)Gfe_i(x), 

where Go = and 

, n 4^ , x q/o , N 2(£;-l) , 

a k (x) = , v (1 - x) 3/2 , 6fc(ar) = 4i (2x - 1) . 

fcv y 2fc + 1 v y w 2fc + l v 7 



(8) 



By induction, 



k / k \ 

Cfc(x) = ^ I ]Q fe s (x) a fe _ i+ i 
*=1 \s=fc-t+2 / 



x) , fc > 1 . 



from which it follows that 

k 



G2(s):=-^L = 4>/2 (!-*)£ 



{ 2 ( 2x-i)}«^} n 



(fc-t)! 1J - 2s + 1 

s=fc-t+l 






r(fc + 3/2) v '^ (fe-t)! 

2 ,/g (fe- 1 )' n T ^y( 2a; )' y r , )t -i- s (*-!)' r(fc-t + 3/2) 



nk + 3/2)" >f^ s\ ^ ' (t-l-s)\ (k-ty. 

o/oJ^lLn ,^(2xr fc ^ 1 (t + g )!r(fc-t- g + i/2) 
= 2 ^f^T3727 (1 - x) §— § ( " 1} -« (k-t-s-iy. • 

It seems impossible to reduce this further. We have 

which combined with (8) gives 

14 7 8 32 

ri = -- + -pi, r 2 = ---p 1 + - P2 . (9) 

The coefficient of pk in the expression for tu is readily seen to be equal to A;!2 fc_1 y / 7r/r(A;+3/2) — > 
2 as k — > oo, by Stirling's formula. 



4 Gram inversion for the stochastic Radon transform and fake 
uniformity 

In this section we consider the stochastic Radon transform, originally introduced in [7] in the 
case where the rotation distribution is Haar measure p, for conjugation- invariant rotation dis- 
tributions, and will utilise results from the previous section. We first briefly repeat the main 
points of the stochastic Radon transform that are relevant for our own purposes; more infor- 
mation can be found in [7] and the references cited therein. The definition of the stochastic 
Radon transform in [7] is based on the projection operator II of an essentially bounded inte- 
grable function p on R 3 , taken to have total mass unity and integral mean zero, with support 
contained in a ball of radius n centred at the origin: 

/•oo 

{n{p}(A))(x 1 ,x 2 ):= Ap(x 1 ,x 2 ,x 3 )dx 3 , A € 50(3) , 



where the rotation group acts on functions on IR 3 via (Ap)(w) = p(A l w). For a random 
rotation, the map II gives rise to a distribution P p on an L 2 space of planar functions 

P p (B) = p{A G SO(3) : U{p}(A) G B} , (10) 

where B ranges over the Borel sets in L 2 space. Note that here we consider what we may call 
the 'population version' of the stochastic Radon transform. 

In [7] it was shown that Pq p = P p for all Q € SO(3), which implies that only the shape of p, that 
is the orbit of p under the above group action, may be recovered in the estimation from sample 
data. This motivates the consideration of the Gram matrix. The Gram matrix of KG M. dxk is 
Gram(y) = V V, where k > 1, and d = 3 for the original shape and d = 2 for its projection. It 
comprises a set of k 'landmarks' as its columns, and is appropriate for describing shape because 
Gvam.{AV) = Gram(V) for any A G SO(3) or SO (2) respectively. 

We now present an extension of the Gram inversion formula from [7, Theorem 4.1, p. 3285] 
to the case of conjugation-invariant, and indeed only shift-symmetric, rotations as defined in 
Section 2. 

Theorem 1. Let V be a 3 x k matrix, k > 1. Let H = diag(l, 1, 0). Then for f P G Z S M , 

f Giam(HAV)f P {A)p(dA) = V T M T (I - D 2 )MV = Gram(F) - Gram(DMK) , 
JsO(3) 

where D 2 = diag ((1 - r 2 )/2, (1 - r 2 )/2, r 2 ). 

Proof. As in [7, pp. 3285-3286], one obtains from Proposition 1 and the fact that R = R that 



■MV 



[ Gram{HAV)f P (A)p(dA) = V T M T ■ \l - E {(.Re (3) )(.Re (3) ) T } 

= Gram(y) - Giam(DMV) , 
D 2 = E{(Re^)(Re^) T } = ([ w t w J f° Re(3) (w)a 2 (dw) 

From Proposition 1 it follows that D 2 is diagonal with E{tr(L> 2 )} = E{(Re^) T (Re^)} = 1. 
This completes the proof. □ 

Remark 2. The first equation of the above proof is easily adaptable to rotations in M rf , d > 2. 

The following example reveals important consequences of the results developed so far. 

Example 1. For the Cayley-LMR distribution defined at (3), from Propositions 1 and 2 we 
rederive for the case of rotations in M 3 the result stated for general finite dimension in [4, 



Proposition 3.3, p. 419]: 



'Re 

Theorem 1 and equations (9), (5) yield 



f£%(w) = 2' K (« + 1) (1 + w 3 ) K , weS 2 . 



[ Gram{H AV}f CAY {MtK) ( A) fi(dA) = Gram(y) - Gram(L>^ AY My) , 
JsO(3) 



ISO{3) 

Z^ AY = diag 



cay - / ; 20c + 1) / 2(k + 1) 12 + k + k 2 



6 + 5k + k 2 V 6 + 5k + k 2 V 6 + 5k + k 



Regarding n, T2 as functions of «, we have T2(k) = 1/3 not only in the Haar measure case k = 0, 
but also for k = 1. For the Fisher-von Mises distribution, we obtain from [4, p. 425] that 



f FvM, a _ 2exp{-2«;} /l - x 



7r{i-o(2Kj — i.i(2Kj) V a; 

where Xo,Xi are modified Bessel functions of the first kind, see Appendix 1 in [5], and from this 
we may again obtain the corresponding values of pi,p2 and then of Ti,T2 from (9). Figure 1 
illustrates the behaviour of T2(k) — 1/3, for either the Cayley-LMR or the Fisher-von Mises 
density, for k G [0, 1]. The instance k = \ for the Cayley-LMR distribution is of special interest 
and motivates the following definition. 

Definition 1. Let V be a family of rotation distributions as described in Section 2. We say 
that V has the faice uniformity property if there exists k > such that the Gram matrix in the 
context of Theorem 1 has the same distribution as for fi. 

The fake uniformity property sets the Cayley-LMR distribution distinctly apart from the Fisher- 
von Mises distribution in a novel aspect, besides the matter of the singular behaviour of the 
former at G = it already noted in [4, p. 425]. It considerably complicates, if not frustrates, 
any attempt at the development of a statistical test for the rotation distribution to be Haar 
measure, if only the random projections are known. For it to be present, a sufficient condition 
is that the derivative of T2(k) at k = 0+ is negative. By (9), the latter condition is satisfied if 
and only if 

^p 2 (0+)<pi(0+). 

Note that this condition would not change if k was instead reparametrised by a monotone 
function k = (p(k) with <f(0) = 0. 

Finally, we note that for both distributions t 2 (k) — > 1, and hence D K — > diag(0, 0, 1) and 



/ Gram(HAV)fp(A)fi(dA) ->• V T ■ \l - E {(Me (3) )(Me (3) ) T ) 
JSO(3) L l J 

Gram {(i- (Me (3) )(Me (3) ) T ) v\ 



■V 



as k — > oo. This is intuitive: as the concentration parameter diverges to infinity, for Gram 
inversion, only the projection to the space orthogonal to the last column of M is relevant. 
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Figure 1: Values of T2(k) — 1/3 as a function of n for the Fisher-von Mises and the Cayley-LMR 
distribution. 



5 Bayes classifier for shift-symmetric rotations 

We continue to use the notation and framework from Section 4. Consider now two shift- 
symmetric random rotations I ^ Pi S Ttf *, i = 1, 2, with Mi 7^ M.%- There exist Q E 0(3) and 
a £ (0, 7r) such that 



I - MxM^ = Q T A{a)Q with A(a) = (o^(a)) 



A 



V 



cos a 



sin a 



(T 



sin a 1 — cos a 
0, 



(11) 



The angle a in (11) is the Riemannian distance between the rotations Pi when coupled as in 
our definitions; see [6]. Our main result in this section is as follows. For simplicity we assume 
equal prior probabilities tt\ = 1x2 = 0.5 for selecting an object distributed as Pi or Pi- 

Theorem 2. With the foregoing notation and conditions, let T denote the event of correctly 
classifying a Class 1 [respectively, arbitrary] observation to Class 1 [respectively, its true class] . 
Then 



P(T) 



p(e<?H*p("£*£* 



+ 



2) 2 

tan (f ) 



2, 



E^cot 



+ 



sin 



e 
2 ~ {""" v 2" 

-r-X E -^ COt — 

(f) I V2 



a 



< 9 < vr 



o 



7T - — < G < 7T 



Remark 3. The Gram matrix is generally unsuitable for classification. To appreciate why, 
take Mi such that M±e^' = M2e^>. The two shapes of the projections do not differ and hence 



10 



classification is impossible. 

Proof of Theorem 2 is at the end of this section. We first discuss classification accuracy for 
small a. Define 

Kx) •■= \hr— fx(x) , o<x<i. 

V 1 — X 

Note that for the Cayley-LMR distribution with M = I, we have h(x) = x k /B(k + 1/2,3/2). 
Regarding the left-hand side in Theorem 2 as a function ipR(a) and writing w a = cos(a/2), we 
have 

V'hC") = / \ h(x)dx + - / J h(x)dx 

l+wa V X 2 l l-w a V X 

2 2 

tan(f) r^P / x 1 r 1 ^ / ^ 

+ -^ I h(x) dx + — - - / h(x) dx . 

It is readily checked that, with the ± signs taken respectively, 

{ tan (l)F 



i- 1 ^ r /a\i±i 



\ 1 -I- 1 ± w o 

V 1 ^ 2 
from which we obtain after some manipulations the simple expression 

— / " /i(x)dx-- / " h(x)dx\ , (12) 

a weighted difference between two integral means. This expression is analogous to the well- 
known one for the case of two multivariate Normal distributions with equal covariance matrix 
given in [10, p. 22], when a = 0. One may also regard it as a measure for deviation from 
uniformity, but unlike the one presented in the previous section it is a non-local one and hence 
probably less useful in practice. 

Proof of Theorem 2. We only need to prove the assertion in the case where T is the event of 
correct classification into Class 1. We have 

P(Observation correctly assigned to Class 1) = P(tr{PM 1 T (I - MiMj)} > 0) . 

Moreover, with the notation at (11), 

tr{PM 1 T (/ - M 1 M%)} = tr {PMjQ T A{a)Q} = ti {Q P Mj Q T A(a)} = ti{PMjA{a)} . 



11 



Using the representation at (2), we have 

PMl = (bij) = (1- cos e)UU T + cos(0) I + sm(0) S(U) , 

b u = (1 - cos Q)Uf + cos 9 , 612 = (1 - cos S)U 1 U 2 - U 3 sin 6 , 

b 21 = (1 - cos 0)U 1 U 2 + U 3 sin 6 , 6 22 = (1 - cos 6)C/ 2 2 + cos 6 . 

We therefore obtain 

2 
tr(PMf^) = Y^ b ijaji(a) = (1 - cosa)(6n + 622) + (sina)(6 2 i - 612) 

= (1 -cosa)(l - cos9)(l - C/f) + 2(1 - cos a) cos 9 + 2(sina)C/ 3 sine 

= aUl + bU 3 + c, 

a = -2(l-cosa)(l-X) <0, 6 = 4(sina)v^(l - *) , c = 2(1 - cosa)X . 

(13) 

(i) / rr (2) .., , r (i,2) 



The zeros of this quadratic in C/3 are C/3 < C/3 with C/3 = (2a) 1 (6 ± V& 2 — 4ac) , where 



sin a I X fa\ f@ 

cot — cot 



2a 1 - cos a V 1 - ^ V 2 / V 2 

Vb 2 - 4ac _ _ v/l^sin 2 a)X(l - X) + 16(1 - cosa) 2 X(l - X) 
2a -4(l-cosa)(l-X) 



sin 2 a + (1 — cos a) 2 cot ( — I = esc ( — I cot I — ) , 



1-cosav <• ' ^2 J V 2 / \ 2 

hence (noting that cot/3 — csc/3 = — tan(/3/2) and cot/3 + esc/? = cot(/3/2)) 

UP = - tan f ^ cot (%\ , C/ q (2) = cot (^ cot / ° 



3 - ^ u;w ^ 2 y , ^3 - U y V 2 

It is easily seen that C/3 < — 1 if and only if B < a/2, and C/3 > 1 if and only if < 7r — a/2. 
Thus the integration with respect to the distribution Fjj 3 of C/3 is over the domains 

1 1 f/ (2) 

«=>- 6 < - , / <S=^ - < < IT , / <=^ 6 > 7T . (14) 

-1 2 JjjW 2 " " 2 JjjW 2 

The density of arcsinC/3 is {7 1— >■ (cos7)/2} for 7 G [— 7r/2, tt/2]. Changing variables and 



12 



integrating over Fjj 3 gives for the three expressions in (14), respectively, P(0 < a/2), 

c-k-ol/2 [-1 i i-w-a/2 f /p.- 

F e (d9) F U3 (du) = - '• ' " 

la/2 Jul 1 ' l Ja/2 



!T Fem L F »> (du > - \ C I 1 + tan ® cot (?) I FaidS) 



L/a „ a\ tan (f If /0\ a ^ a 

_ P (_<e< 7r --) + ^i2 E | cot ^_j ; _<e< 7 r-- 
^ 2 > ,_ i r /e\ 



F @ (d9) / F [/3 (d«) = ^ T ^ Y E cot - ;7r--<e<7r , 



where we used the fact that 



S 1) = {tan(f)+cot(f)}cot(f) 



^ - ^ - {tan («) + cot ( «) | cot ( ? ) = ^ cot ( @ 



(f) V2 
This concludes the proof of Theorem 2. □ 
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